Reverse-time migration-based reflection tomography using teleseismic free surface multiples

Converted and multiply reflected phases from teleseismic events are routinely used to create structural images of the crust-mantle boundary (Moho) and the elasticity contrasts within the crust and upper mantle. The accuracy of these images is to a large extent determined by the background velocity model used to propagate these phases to depth. In order to improve estimates of 3-D velocity variations and, hence, improve imaging, we develop a method of reverse-time migration-based reflection tomography for use with wavefields from teleseismic earthquakes recorded at broad-band seismograph arrays. Reflection tomography makes use of data redundancy – that is, the ability to generate numerous structural images of the subsurface with different parts of the wavefield. In exploration seismology (where it is known as migration velocity analysis) reflection tomography typically involves the generation of an extended image (e.g. offsetor angle-gathers), and the fitness of the background model is evaluated through the application of image-domain annihilators. In regional-scale passive source seismology, however, annihilation-based methods are inadequate because the sparse and irregular distribution of teleseismic sources is not likely to produce illumination over a sufficient range of angles. To overcome this problem we turn towards a source-indexed moveout scheme. Instead of extended image annihilation, we determine the success of the tomographic velocity model by cross correlating images produced with multiply-scattered waves from different teleseismic sources. The optimal velocity model is the one that minimizes correlation power between windowed images away from zero depth shift. We base our inversion scheme on the seismic adjoint method and a conjugate gradient solver. For each image pair, the update direction is determined by correlations between downgoing wavefields with upgoing adjoint wavefields for both images. The sensitivity kernels used in this method is similar to those found in other forms of adjoint tomography, but their shapes are controlled by the spatial distribution of the error function. We present the method and a proof-of-concept with 2-D synthetic data.